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ABSTRACT 



We present the analysis of the first f 8 months of data obtained with the COS- 
MOSOMAS experiment at the Teide Observatory (Tenerife). Three maps have been 
■ obtained at 12.7, 14.7 and 16.3 GHz covering 9000 square degrees each with a reso- 

O-f lution of ~ 1 degree and with sensitivities 49, 59 and 115 /iK per beam respectively. 

q • These data in conjuction with the WMAP first year maps have revealed that the 

?H ' Cosmic Microwave Background (CMB) is the dominant astronomical signal at high 

^ \ galatic latitude (|6| > 40°) in the three COSMOSOMAS channels with an average 

d • amplitude of 29.7 ± 1.0 (68% c.l. not including calibration errors). This value 

is in agreement with the predicted CMB signal in the COSMOSOMAS maps using 
the best fit A-CDM model to the WMAP power spectrum. Cross correlation analysis 
of the 408 MHz map and the COSMOSOMAS data at high galactic latitudes give 
J-j | values in the range 17.0 — 14.4 fj,K from 12.7 to 16.3 GHz. Removing detected point 

sources in this template, reduces the amplitude of the correlated signal to 8-9 ^tK. The 
mean spectral index of the correlated signal between the 408 MHz desourced and the 
COSMOSOMAS maps is between -3.20 and -2.94 at |6|>40° which indicates that this 
signal is due to synchrotron emission. Cross-correlation of COSMOSOMAS data with 
the DIRBE map at 100 /im shows the existence of a common signal with amplitude 
7.4 ± 1.1, 7.5 ± 1.1, and 6.5 ± 2.3 /*K in the 12.7, 14.7 and 16.3 GHz COSMOSO- 
MAS maps at |6| >30°. Using the WMAP data we find this DIRBE correlated signal 
rises from high to low frequencies flattening below ~ 20 GHz. At higher galactic lati- 
tudes the average amplitude of the correlated signal with the DIRBE maps decreases 
slightly. The frequency behaviour of the COSMOSOMAS/WMAP correlated signal 
with DIRBE is not compatible with the expected tendency for thermal dust. A study 
of the Ha emission maps do not support free-free as a major contributor to that signal. 
Our results provide evidence of a new galactic foreground with properties compatible 
with those predicted by the spinning dust models. 

Key words: cosmology: observations - cosmic microwave background - galactic 
anomalous emission 



1 INTRODUCTION 

In recent years many experiments have measured Cos- 
mic Microwave Background (CMB) anisotropics at angu- 



lar scales from several degrees to a few arcmin. The anal- 
ysis of the resulting power spectra has tested the Lambda 
Cold Dark Matter (ACDM) model and allowed an accu- 
rate determination of cosmological parameters (e. g. Spergel 
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et al. 2003; Rebolo et al. 2004). Detailed knowledge of the 
galactic foregrounds (synchrotron, free-free and vibrational 
thermal dust emission) is required to ensure the validity 
of the results. The joint analysis of different templates of 
galatic emission with WMAP observations (Bennett et al. 
2003; Hinshaw et al. 2003) has allowed a better assesment 
of the properties of these galactic foregrounds at the frequen- 
cies where CMB searches are conducted (see also Banday et 
al. 2003). Apart from the three classical diffuse components 
(synchrotron, free-free and vibrational dust emission), there 
is increasing evidence of an additional dust correlated fore- 
ground (Casasus et al. 2004; de Oliveira-Costa et al. 1997, 
1998, 1999, 2004; Finkbeiner et al. 2002; Kogut et al. 1996) 
which does not follow the expected thermal dust spectrum 
at frequencies below v <40 GHz. The physical nature of such 
a component is intriguing (Lagache 2003; Leitch et al. 1997; 
Lim et al. 1996; Mukherjee et al. 2001; Banday et al. 2003). 
Although it was initially ascribed to free-free emission, the 
low level of related Ha emission makes unlikely this expla- 
nation. Alternatively, one of the most attractive scenarios 
to explain the origin of that emission is electric dipole ra- 
diation from spinning small dust particles as proposed by 
Draine & Lazarian (1998a, b). The key observational pre- 
diction of such a model is a turn-over in the spectrum at 
frequencies 10-20 GHz which are not covered by COBE or 
WMAP. 

The Tenerife experiment (Gutierrez et al. 2000) has pro- 
vided observations at 10 and 15 GHz able to test this pre- 
diction (e.g. de Oliveira et al 2004). The COSMOSOMAS 
(COSMOlogical Structures On Medium Angular Scales) ex- 
periment represents a qualitative improvement in terms of 
angular resolution, sensitivity and sky coverage to study any 
possible new foreground in the frequency range 10-17 GHz. 
This experiment, designed and built at the Instituto de As- 
trofi'sica de Canarias, is located at an altitude of 2400 m 
in the Teide Observatory, Tenerife (Spain). It consists of 
two similar instruments working at frequencies centred at 
11 and 15 GHz named COSMO 11 and COSMO 15 respec- 
tively. This paper is dedicated exclusively to the analysis 
of COSMO 15 data. The COSMO 11 results will be pre- 
sented in a separate paper (Hildebrandt et al. 2005 in prep.). 
A detailed technical description of the COSMO 15 instru- 
ment has been presented in Gallegos et al. (2001) where we 
also reported the first results obtained with this experiment. 
Here, we present additional data taken up to January 2001 
improving the sensitivity by a factor two. We use the COS- 
MOSOMAS data in conjuction with the first year WMAP 
data and several Galactic templates to make a robust es- 
timation of the CMB and the different components of the 
galactic diffuse emission in the frequency range 12-17 GHz. 



2 EXPERIMENTAL SETUP, OBSERVATIONS 
AND DATA REDUCTION 

2.1 Instrument 

COSMO 15 uses a flat rotating circular mirror (2.5 m of di- 
ameter) made of aluminum canted by 5 degrees with respect 
to its rotation axis. The sky radiation is reflected by that 
mirror and redirected to a 2.4 m off axis parabolic antenna 
where it is focused and detected by a cryo-cooled HEMT am- 
plifier. The rotation of the mirror produces circular scans in 



Table 1. The three channels of the COSMO 15 instrument. 



Channel 


Freq. 


Bandwidth 


T sys 


FWHM HA 


FWHM Dec . 




(GHz) 


(GHz) 


(K) 


(deg) 


(deg) 


CI 


12.7 


1.17 


38.5 


0.92 ± 0.01 


1.11 ± 0.01 


C2 


14.7 


0.94 


33.5 


0.81 ± 0.01 


1.00 ±0.01 


C3 


16.3 


1.18 


37.9 


0.79 ± 0.03 


0.85 ± 0.02 



Table 2. Observations analysed in this work 



Epoch 


Range in Decl. 


CI 


C2 


C3 


(dd/mm/yy) 


(deg) 


(Number of days) 


27/09/99-08/11/99 


16.7-36.7 


3 


4 


5 


01/12/00-15/02/00 


23.7-43.7 


13 


10 


7 


20/02/00-15/03/00 


27.4-47.4 


10 


9 


7 


22/03/00-18/09/00 


24.8-44.8 


56 


56 


45 


06/11/00-16/01/00 


26.4-46.4 


36 


27 


20 



the sky at 1 Hz. The data are sampled ~ 220 times per turn 
of the mirror. Earth rotation allows the daily observation of 
a sky strip complete in right ascension and 20 degrees width 
in declination. The instrument has three independent chan- 
nels (hereafter called CI, C2 and C3) working at frequencies 
centred at 12.7, 14.7 and 16.3 GHz respectively with band- 
widths ~ 1 GHz and angular resolution ~ 1 degree. The 
instrument is shielded to protect against spillover. To cor- 
rect for changes in the gain of the instrument, every 30 s 
a constant signal at 2 K is injected. System temperatures 
are ~ 35 K with small differences between channels. Ta- 
ble Q summarizes the main parameters including beam sizes 
estimated as explained in Section 3.2. 

2.2 Observations and data procesing 

The data analysed here correspond to the period September 
1999- January 2001. The experiment was running continu- 
ously except for periods of testing, maintenance and tech- 
nical failures, or when the atmospheric conditions were ex- 
treme (strong winds or rain). After careful inspection of the 
data, we selected ~ 100 days of good-quality for the analyses 
of this paper. Table [5] presents a summary of the observa- 
tions indicating the number of useful maps selected in each 
one of the five different instrumental settings of COSM015 
in this period. 

The basic steps in data processing include filtering out 
data affected by technical failures or proximity of the Sun 
and Moon, removal of differential atmospheric emission, ab- 
solute calibration and construction of a daily map for each 
of the three channels. The first steps of data processing are 
done in the Fourier space as explained in Gallegos et al. The 
dominant component in each scan is a spin-synchronous sig- 
nal due to a combination of differential atmospheric emission 
and ground pick-up. As shown in Gallegos et al. (2001), this 
component appears in the data as large scale modulations, so 
it can be accurately removed using a combination of low or- 
der sinusoidal functions. In particular, for the map-making 
procedure of this paper, we have subtracted 7 harmonics. 
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RA and 20 degrees width in declination is obtained for each 
channel. 

The calibration process implies the determination of the 
angular response and the conversion from counts to physical 
units of temperature. A complete calibration of the exper- 
iment is performed on a daily basis using the same stan- 
dard astronomical sources used in the pointing determina- 
tion; Cyg A and Tau A. Both are stable and sufficiently 
strong to be used as absolute calibrators of the daily maps. 
Their fluxes at the COSMOSOMAS frequencies have been 
computed following the model by Baars et al. (1977). This 
model has been checked against the fluxes measured for 
these sources by WMAP 1 in the K, Ka, Q and V chan- 
nels, and with the direct measurement of Cyg A at 15 GHz 
by Baker et al. (1977). We estimate the uncertainty in cal- 
ibration of our data is ^ 10 %. The FWHM of the beam is 
determined from the radial profile in Dec. and RA. of these 
calibration sources. Average values of the FWHMs measured 
in all daily maps and their errors are listed in Table Q The 
conversion factor is then established from source flux and 
beam FWHM and apparent amplitude in counts. In the cal- 
ibrated daily maps we find typical sensitivities are 0.5-1.0 
mK per beam. 

3.2 Noise analysis 

A robust estimation of the nature and amplitude of the as- 
tronomical signals present in the COSMOSOMAS maps re- 
quire an accurate knowledge of the statistical properties of 
the noise. The noise in each daily map is characterised by the 
noise covariance matrix between each pair of pixels. A rea- 
sonable first approximation assumes that this only depends 
on angular distance. This is further simplified if the noise is 
white, in which case the covariance matrix is diagonal. In our 
case, the noise matrix can be estimated from the analysis of 
data at high galactic latitude in single day observations. In 
such daily maps the relative contribution of the astronomical 
signals is small compared to the dominant noise component. 
Figure shows the correlation functions (assuming that the 
covariance between any pair of pixels is only a function of 
the distance) evaluated from one typical day of data in the 
three channels of COSMOSOMAS. For this analysis, we se- 
lected the region enclosed in a square RA=(115°, 267°) and 
Dec. =(25°, 45°), excluding also the pixels affected by the 
presence of radiosources according to the mask described 
in Section 4. The presence of non diagonal elements in the 
covariance matrix is revealed by non-zero values in the cor- 
relation function at angles different from zero. Whether this 
is due to the instrument or to the atmosphere is analysed 
below. 

Because the observations are simultaneous in the three 
channels, some degree of correlated signal between them is 
expected due to common parts in the receiver system and 
contamination by atmospheric signals. We have estimated 
this by cross correlating the simultaneous observations in 
the three channels. Table [3] presents the average noise val- 
ues (diagonal) and square root of the average covariances 

1 The WMAP analysed in this paper 

are the first year full sky maps (see 
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Figure 1. Effective window function of the COSMOSOMAS 
maps. 



Further additional components due to changes in the atmo- 
spheric conditions during the day are also eliminated by a 
long period fit to each point of the scans over the day. These 
two processes completely eliminate the response of the ex- 
periment to temperature fluctuations on scales ^ 5 degrees. 

We follow a numerical approach to evaluate the re- 
sponse of the experiment to different angular scales. This is 
done simulating the reponse to a point source after passing 
through the scanning strategy of the experiment. Given that 
the data have been obtained with five different instrumental 
settings (see Table |5J, we have defined an effective window 
function as a weighted mean value of the window functions 
associated with each setting. For the different window func- 
tions the integral under the function is used and for different 
pointings the weight is proportional to the number of days 
observed. Figure Q presents this effective window function 
of the experiment where we have taken into account that 
the final maps are convolved to a common resolution of 1.12 
degrees (see next sections). 



3 MAPS 

3.1 Daily maps and astronomical calibration 

Using the strong astronomical radio sources Cyg A 
(RA=299.87°, Dec.=40.73°, Eq. 2000) and Tau A 
(RA=83.64°, Dec. =22. 01°, Eq. 2000) we accurately obtain 
the optical configuration which in turn is used to determine 
the path swept out by the circular scanning. With this the 
celestial coordinates of data in each position of the scans 
can by calculated and used to project them into a grid with 
| x | degrees in RA and Dec. This grid was choosen to ade- 
quately sample the beam. A daily map with full coverage in 



'://lambda.gsfc. nasa.gov/product/map/drl/imaps_frequency.cfm 
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Table 3. Average values for the noise and cross-correlations per 
pixel in the daily COSMOSOMAS maps. 





CI 


C2 


C3 




(mK) 


(mK) 


(mK) 


CI 


2.11±0.73 


0.50 ±0.10 


0.42 ±0.13 


C2 




2.03±0.92 


0.42 ±0.11 


C3 






3.30±1.12 



between daily maps in each channel over the whole period 
of observations. The listed values are derived by the follow- 
ing expression 

rpk rpk' 

Cov(k,k') = ±^±^ gJ-Zl , k,k' = C1,G2,G3 

where Ty corresponds to the temperature of the band k (= 
CI, C2, C3), in pixel i and in map j (j runs over all daily 
maps). This correlated signal between channels extends up 
to angular scales of 5 degrees. We attribute this daily correla- 
tion to the presence of atmospheric noise which still persists 
after the filtering of the low order harmonics. 

We have checked for the presence of possible correla- 
tions between consecutive daily maps and conclude that the 
correlation is below 1 %. This is at the level expected from 
astronomical signals. We conclude that the atmospheric 
noise in the individual COSMOSOMAS maps is uncorre- 
cted from day to day, and thus each daily map is considered 
as independent in the stacking of the data. 

3.3 Stacked maps 

For each channel, the n daily maps are stacked together 
to get final maps with improved signal to noise ratio by a 
factor y/n. These final maps are presented in Figure They 
are complete in right ascension and cover from 16.7 to 47.4 
degrees in declination. The conspicuous structures at RA~ 
300 degrees correspond to the main crossing of the galactic 
plane. The two strong sources Tau A and Cyg A are clearly 
visible. Other sources are also present in the maps and will 
be described below. The regions with higher sensitivity are 
concentrated in the band 25° ^ Dec. ^ 45°. At Dec. ^ 
16.7° the number of observations is smaller and the noise 
is higher as is easily appreciated in Figure [3] The errors 
are calculated from the dispersion of all the measurements 
contributing to a given pixel of the map. The mean noise 
in a wide region (120° < RA < 270°, 25° ^ Dec. < 45°) 
excluding the galactic plane, are in 173, 185 and 329 ^iK 
per pixel for channels CI, C2 and C3 respectively; these are 
translated into 49, 59 and 115 fiK beam -1 . The rms of the 
three maps computed in the same region are 90, 90 and 116 
/j,K beam - . These larger values of the rms as compared with 
the noise indicates the presence at high galactic latitudes 
of clear statistical signals in each map. The evaluation and 
origin of these signals is considered in the following sections. 

Figures 0| and show a visual comparison between the 
COSMOSOMAS maps at the three observing frequencies, 
and the WMAP map at 23 GHz (all convolved to a common 
circular resolution of 1.12°, which corresponds to the worse 
resolution of those presented in Table 1). The two regions 
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Figure 2. Correlation function of one day of COSMOSOMAS 
data at the three observing frequencies. 



displayed are located in the galactic plane, and the promi- 
nent structures are mainly due to diffuse emission from our 
galaxy. The remarkable agreement between the structures 
observed by COSMOSOMAS and WMAP supports the re- 
liability of the former. Table 0] lists the observed fluxes of 
some of the strongest radio sources identified in the COS- 
MOSOMAS maps. The fluxes are determined by the same 
method used in Watson et al (2005) where the best fitting 
amplitude for the expected modelled source profile to the 
data is found. For a comparison Table [I] presents also the 
fluxes measured for these sources in the two low frequency 
WMAP maps (Bennett et al. 2003). COSMOSOMAS and 
WMAP data are not coeval, and so this comparison is only 
indicative for variable sources (3C84, 4C39.25, 3C345, etc). 
As expected in the frequency range covered by COSMOSO- 
MAS and WMAP most of the sources have a flat spectrum. 
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Table 4. Fluxes of several point sources identified in the COSMOSOMAS CI, C2 and C3 maps. WMAP fluxes in the K and Ka channels 
(taken from Bennett et al. 2003) are also listed. Errors and upper limits are 68 % and 95 % c.l. respectively. 



COSMOSOMAS WMAP 
ID RA, Dec F C i F C2 F C3 F K F Ka 





(J2000) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


0234+28 


(39.5°, 28.8°) 


2.6 + 0.3 


4.1 + 0.5 


6.0 + 1.0 


4.6 + 0.5 


5.0 + 0.6 


0316+41 (3C84) 


(50.0°, 41.5°) 


17.4 + 1.8 


19.6 + 2.0 


16.8 + 1.9 


12.3 ± 1.2 


8.3 + 0.9 


0433+29 


(69.3°, 29.7°) 


7.7 + 0.8 


8.9 + 1.0 


9.8 + 1.0 


5.3 + 0.5 


4.3 + 0.5 


0738+31 


(115.3°, 31.2°) 


2.4 + 0.4 


1.0 + 0.4 


< 1.4 


1.6 + 0.2 


1.0 + 0.3 


0923+39 (4C39.25) 


(141.8°, 39.0°) 


9.1 ± 1.0 


9.5 ± 1.0 


7.0 + 1.2 


7.8 + 0.8 


6.4 + 0.7 


1328+30 (3C286) 


(202.8°, 30.5°) 


2.4 + 0.3 


3.0 + 0.4 


2.7 + 0.9 


1.5 + 0.1 


0.7 + 0.2 


1640+33* 


(250.0°, 33.5°) 


1.2 + 0.3 


1.6 + 0.4 


1.9 + 0.9 


1.6 + 0.2 


3.1 + 0.3 


1641+39 (3C345) 


(250.7°, 39.8°) 


13.1 ± 1.3 


14.0 + 1.4 


12.8 + 1.5 


8.2 + 0.8 


6.9 + 0.7 


1611+343* 


(243.4°, 34.2°) 


4.4 + 0.5 


4.5 + 0.6 


4.2 + 0.9 


4.5 + 0.5 


4.7 + 0.5 




I 

I 




direct separation of the various astronomical components 
present in the data, namely CMB and the different galactic 
foregrounds. To evaluate the presence of each astronomical 
signal we have applied cross-correlation techniques between 
the COSMOSOMAS maps and existing templates consid- 
ered as good tracers of the different galactic foregrounds. 
These cross-correlation techniques allow an estimation of 
which fraction of the total signal contained in a given map 
is correlated with a reference template. The main assump- 
tion is that the spatial structure in the reference and the 
map under analysis is the same apart from a scaling factor 
in amplitude. 
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Figure 3. Stacked COSMOSOMAS maps at 12.7 (top), 14.7 
(middle) and 16.3 (bottom) GHz. 



The contributions of unresolved sources will be quantified in 
Section 4. 



4 CROSS CORRELATION ANALYSIS 

The limited frequency range spanned by the three COSMO- 
SOMAS channels analysed here, do not allow an accurate 



4.1 Templates 

The maps used as templates are the 408 MHz (Haslam et 
al. 1982), 1420 MHz (Reich & Reich 1986), DIRBE 2 , and 
Ha (Finkbeiner 2003). At 408 MHz the dominant diffuse 
foreground is synchrotron, however the Haslam et al. map 
at this frequency is severely contaminated by the contribu- 
tion of point sources, so as a best tracer of the synchrotron 
emission we have used the 'cleaned map' (desourced) pro- 
vided on the NASA LAMBDA site 3 . The main properties of 
these templates are summarized in Table [K] The 1420 MHz 
has a contribution of both synchrotron and free-free emis- 
sion. Ha emission is a tracer of the warm ionized medium 
which also radiates via free-free emission, although an accu- 
rate estimation of this foreground requires a correction for 
the extinction by dust and the electron temperature of the 
medium (Dickinson et al 2004). However, as we are inter- 
ested in correlated signals at high galactic latitudes where 
the dust extinction is negligible, we have not attempted this 
correction. We only use the DIRBE maps at 100 and 240 
/im as good tracers of galactic dust. We include in the cross- 
correlation analysis also the first year WMAP maps, which 
are dominated by CMB signal, although depending on the 
frequency considered different degrees of galactic contami- 
nation exists. 

Before the correlation analyses, the WMAP and galactic 

2 http:/ /lambda. gsfc.nasa.gov/product/cobe/ 

3 http: / /lambda. gsfc.nasa.gov/product/foreground/haslam_408.cfm 
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Figure 4. From top to bottom: the Galactic plane region ob- 
served by the three COSMOSOMAS channels, and the lowest 
frequency channel of WMAP (before and after filtering according 
to the observing strategy of COSMOSOMAS). 
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Figure 5. As in the previous figure, but for a field towards the 
Auriga-Perseus region. 



template maps have been filtered to the common COSMO- 
SOMAS resolution (1.12°), and have been convolved with 
the simulated circular lock-in analysis and pixelated in the 
same 1/3° grid used before. In the construction of these tem- 
plates small differences will arise with respect to the COS- 
MOSOMAS maps due the 3-sigma clipping algorithm im- 
plemented to reduce the effect of very strong radio sources 
in the observed scans. These differences will be larger as we 
get closer to the galactic plane so we will consider here only 
correlation results at high galactic latitudes \b\ > 30°. Due 
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to the comparatively high amplitude of the galactic plane, 
we do not discard the possibility of some residual systematic 
effects in the band 30° < |6| < 40°. 

4.2 Method 

The cross-correlation between two maps is derived as fol- 
lows. Let d be the data vector containing TV data points. 
We want to fit simultaneously n t "template maps". Let ti 
be the templates (i = 1, ...,n t ), and a a vector with n t ele- 
ments containing the cross-correlation variables in which we 
are interested. Then, the a vector can be derived by mini- 
mizing 

X 2 = (d-a ■ M) T ■ CT 1 ■ (d-a ■ M) (f) 

where C is the covariance matrix including both the signal 
and the noise, and M is a n t x N matrix built from the 
template maps, and whose elements are defined as Mi t j = 
(ti)j, where i = l,...,n t and j = 1, N. Solving for a we 
obtain: 

a = {dC~ 1 M T ){MC~ 1 M T y 1 (2) 

and the corresponding covariance matrix of the a. pa- 
rameters is 

Cov(a) = (MC' 1 M T )~ 1 (3) 

The variance for each parameter ocj is given by the square- 
root of the corresponding element in the diagonal of the 
covariance matrix. Gorski et al. (f996) described this same 
method but in harmonic space. Note that for the case n t = 
1 we recover the standard equations presented in, e.g. de 
Oliveira-Costa et al. (f999), for the case of a single template. 

We have conducted several tests to study the correlated 
component of the noise and the CMB in the covariance ma- 
trix. Our conclusion is that both terms are small as com- 
pared with the white component of the noise and with the 
dominant foreground of each template, and so they can be 
ignored except when WMAP maps are considered as tem- 
plates, in which case it is necessary to take into account the 
CMB component. The analysis has been performed in dec- 
lination strip 25°-44.7° for galactic latitude cuts \b\ > 30, 
40 and 50 degrees to study the dependence of the different 
contributors with galactic latitude. 

We note that the harmonic supression procedure fol- 
lowed in our processing of the COSMOSOMAS data re- 
moved a large proportion of the correlations induced by the 
large scale structure of our galaxy (e.g. cosec law distribu- 
tion). 

We assume a model in which each COSMOSOMAS and 
WMAP map contains a superposition of CMB, diffuse galac- 
tic fluctuations, and noise. Due to the high sensitivity of the 
WMAP maps, their structures are dominated largely by as- 
tronomical signal. The WMAP collaboration provides sev- 
eral masks which account for sky regions near the Galactic 
plane or intense radio-sources. We have used the so-called 
KpO mask (Bennett et al. 2003) which eliminates 2f % of the 
area of the sky, mostly situated in the galactic plane. This 
mask excludes also a region within 0.6° around the ~ 700 



Table 5. Galactic templates The third column indicates the dom- 
inant components in each template (1. CMB; 2. Synchrotron; 3. 
Free-free; 4. Dust). 



Map 


Frequency 


Component 


Sky coverage 




Haslam et al. (1982) 


408 MHz 


2 


Full sky 


Reich &: Reich (1986) 


1420 MHz 


2, 3 


-19° ^ 5 ^ +90° 


DIRBE 


100, 240 ftm 


4 


Full sky 


Finkbeiner (2003) 


H a 


3 


Full sky 


WMAP 


23, 33, 41, 61, 94 GHz 


1, 2, 3, 4 


Full sky 




RA (degrees) 



Figure 6. Mask of discrete sources and galactic emission used in 
the correlation analysis. 



point sources contained in the catalogue by Stickel et al. 
(f994), and in catalogues of AGNs and X-ray sources (Ben- 
nett et al. 2003). We have also masked an additional source 
located at RA=250.33°, Dec.=33.66° (see Table 4) which 
has ~ 1.4 and ~ 0.5 Jy in the CI and in WMAP_K maps 
respectively. To account for the large beam of COSMOSO- 
MAS as compared with WMAP, we have extended the KpO 
mask by 1 pixel (0.33° x 0.33°) surrounding all of the masked 
sources. For the three most intense sources (3C45, 4C39.25 
and DA406) the mask was extended to cover a radius of 
3 degres, in order to remove scanning artifacts larger than 
the noise. Figure HJ presents this mask. 

4.3 Results 

We have obtained all the correlations by fitting a single tem- 
plate, except for the case of cross-correlations with infrared 
and Ha maps where both maps have been fitted simultane- 
ously since we expect a significant coupling between them. 
Indeed, for the case of \b\ > 30°, we find a coupling of 18 % 
for Ha-DirbelOO /im, and 8 % for Ha-Dirbe240 /im. 

The cross-correlation results between WMAP and COS- 
MOSOMAS are shown in Figure □ and Table © while the 
results of the correlations with the 408 MHz, 1420 MHz and 
Ha templates are summarized in Table |7| and Figure [S] The 
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Table 6. Correlated signal between the WMAP and COSMO- 
SOMAS maps for galactic latitude cuts (from top to bottom) of 
[6| > 30°, > 40° and > 50°. Numbers correspond to the product 
of the coupling constant a and the rms of the template WMAP 
map ctwmap- All values are in /zK. 

WMAP COSMOSOMAS 



Channel 


a WMAP 




Cl 




C2 




C3 








>> 


> 30° 












WMAP.K 


33.6 




35.4 ± 1 


1 


34.3 ± 1 


2 


38.3 ± 2 


5 


WMAP.Ka 


30.1 




31.5 ± 1 


1 


28.4 ± 1 


2 


31.9 ± 2 


5 


WMAP.Q 


30.0 




29.7 ± 1 


2 


27.6 ± 1 


2 


31.0 ± 2 


5 


WMAP.V 


29.0 




27.3 ± 1 


1 


25.9 ± 1 


2 


26.3 ± 2 


5 


WMAP.W 


28.6 




26.3 ± 1 


1 


25.3 ± 1 


2 


27.4 ± 2 


5 






>> 


> 40° 












WMAP.K 


33.5 




38.6 ± 1 


4 


36.3 ± 1 


5 


42.5 ± 2 


9 


WMAP.Ka 


30.0 




34.1 ± 1 


4 


30.0 ± 1 




35.9 ± 2 


9 


WMAP.Q 


30.0 




32.4 ± 1 


4 


30.0 ± 1 


5 


34.8 ± 2 


9 


WMAP.V 


29.0 




30.7 ± 1 


4 


28.2 ± 1 




30.2 ± 2 


9 


WMAP.W 


28.7 




29.8 ± 1 


4 


27.9 ± 1 


5 


32.4 ± 2 


9 






M 


> 50° 












WMAP.K 


33.6 




38.0 ± 1 


5 


33.8 ± 1 


6 


42.3 ± 3 


2 


WMAP.Ka 


30.3 




35.3 ± 1 


5 


28.4 ± 1 


6 


36.0 ± 3 


2 


WMAP.Q 


30.6 




33.1 ± 1 


5 


28.8 ± 1 


6 


36.3 ± 3 


2 


WMAP.V 


29.5 




31.4 ± 1 


5 


26.8 ± 1 


6 


31.1 ± 3 


2 


WMAP.W 


29.0 




31.5 ± 1 


5 


28.0 ± 1 


6 


33.7 ± 3 


2 




40 60 
Frequency (GHz) 



i 



|b| > 30" 
|b| > 40" 
|b| > 50' 



6 b 

|b| > 30' 
|b| > 40" 
|b| > 50" 



80 

|b| > 30* 
|b| > 40" 
|b| > 50' 



40 



80 



Frequency (GHz) 



Figure 7. Common signal between COSMOSOMAS (top : 12.7, 
middle : 14.7, bottom : 16.3 GHz) and WMAP for three galactic 
cuts. 



results of the cross-correlation analysis with the dust tem- 
plates (100 fim and 240 /im) are summarized in Table [8] 
There are significant correlations between the three chan- 
nels of COSMOSOMAS and the two DIRBE frequencies. 
Figure[5]shows the signal correlated with the map at 100 fim. 
The results quoted in the previous tables and figures include 
only the uncertainty due to statistical errors. As a realistic 
test on the significance of these errors, we performed cor- 



Figure 8. COSMOSOMAS and WMAP correlated with the syn- 
chrotron maps at 408 (top) and 1420 MHz (bottom) for three 
galactic cuts. 



relations of the different templates with randomly ordered 
pixel values and maps rotated with respect the Galactic axis. 
While the random pixel ordenations basically confirm the 
statistical errors, the dispersion of the correlations of the 
rotated maps were typically 2 times larger. Thus, we will 
consider as significant those correlations with amplitude at 
least 2 times the statistical error. 



5 DISCUSSION 

We find a highly significant correlated signal between the 
COSMOSOMAS and WMAP maps in the three galactic 
cuts. The amplitude of the common signal ranges between 
~ 30 and ~ 40 //K depending of the pair of channels and 
the galactic cut considered. The amplitude of this signal 
decreases as we move from low to high frequency WMAP 
channels. This is a clear indication of the presence at some 
level of foreground signals with a negative spectral index. 

The average values found in the correlation between the 
COSMOSOMAS maps and the WMAP channel is 29. 7± 1.0 
/xK, which corresponds to the weighted-mean of the values 
found between WMAP_V and the three frequencies of COS- 
MOSOMAS for a galactic cut |6| > 40° (the V band was 
choosen because it presents the lowest galactic contamina- 
tion). This value is in good agreement with the amplitude 
of the signal in the filtered WMAP.V channel (29.0 at 
\b\ > 40°), and with the amplitude of the CMB signal ex- 
pected in the COSMOSOMAS maps according to simula- 
tions using the best fit A-CDM model to the WMAP power 
spectrum (« 30 jjK). This result indicates the realibility of 
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Table 7. Correlated signal between the COSMOSOMAS - WMAP maps and three galactic templates. All correlation values have units 
of fiK, and were derived as the product of the coupling constant a times the rms of the corresponding template (rrxemplate)- Values for 
408 MHz and 1420 MHz maps were obtained using a single-template fitting, while the values for H« were obtained from a simultaneous 
fitting with the Dirbe lOOfim map. For the 408 MHz map, we also consider the "desourced version" (labelled as Dss, see text for details). 



Template 




° 'Tempi a 


ta 


CI 


C2 


C3 


WMAP.K 


WMAP.Ka 


WMAP.Q 


WMAP.V 


WMAP.W 




\b\> 30° 


408 MHz 


4 


88 X 10 5 


fj,K 


17.0 ±1.1 


12.3 ± 1.2 


15.4 ± 2.4 


4.7 ± 0.3 


2 


1 ± 0.3 


1.7 ± 


.3 


1.1 ± 0.3 


0.6 ± 0.3 


408 MHz (Dss) 


4 


79 X 10 5 




9.3 ± 1.1 


8.7 ± 1.2 


7.3 ± 2.4 


3.7 ± 0.3 


2 


± 0.3 


1.8 ± 


.3 


1.4 ± 0.3 


1.1 ± 0.3 


1420 MHz 


2 


54 X 10 4 


IJ.K 


20.7 ±1.1 


13.7 ± 1.2 


13.2 ± 2.5 


5.2 ± 0.3 


2 


1 ± 0.3 


1.4 ± 


.3 


0.6 ± 0.3 


0.0 ± 0.3 


Ha 




0.07 R 




2.6 ± 1.1 


1.4 ± 1.2 


-2.2 ± 2.4 


0.1 ± 0.2 





5 ± 0.2 


0.1 ± 


.2 


0.1 ± 0.2 


0.4 ± 0.2 


\b\> 40° 


408 MHz 


4 


74 X 10 5 


IJ.K 


17.0 ± 1.3 


13.4 ±1.4 


14.4 ± 2.8 


6.1 ± 0.3 


2 


8 ± 0.3 


2.3 ± 


.3 


1.6 ± 0.3 


1.1 ± 0.3 


408 MHz (Dss) 


4 


71 X 10 5 




8.0 ±1.4 


8.9 ±1.4 


9.4 ± 2.8 


3.S ± 0.3 


2 


2 ± 0.3 


1.6 ± 


3 


1.4 ± 0.3 


0.9 ± 0.3 


1420 MHz 


2 


48 X 10 4 


fj.K 


19.1 ± 1.4 


13.4 ±1.5 


12.1 ± 2.9 


8.1 ± 0.2 


2 


7 ± 0.3 


1.9 ± 


.3 


1.2 ± 0.3 


0.6 ± 0.3 


Ha 




0.06 R 




1.1 ± 1.4 


5.1 ± 1.5 


1.0 ± 2.9 


0.3 ± 0.2 


(") 


8 ± 0.2 


0.0 ± 


. 2 


0.2 ± 0.2 


0.4 ± 0.2 


\b\> 50° 


408 MHz 


4 


59 x 10 5 


li K 


13.8 ± 1.5 


9.4 ± 1.6 


10.5 ±3.1 


5.8 ± 0.3 


3 


5 ± 0.3 


2.8 ± 


.3 


1.9 ± 0.3 


1.7 ± 0.4 


408 MHz (Dss) 


4 


69 X 10 5 


IJ.K 


5.5 ±1.5 


5.6 ±1.6 


6.0 ±3.1 


3.4 ± 0.3 


2 


2 ± 0.3 


1.7 ± 


3 


1.2 ± 0.3 


1.0 ± 0.4 


1420 MHz 


2 


45 x 10 4 




15.3 ± 1.6 


9.6 ± 1.6 


6.8 ± 3.3 


5.8 ± 0.3 


3.0 ± 0.3 


2.0 ± 


.3 


1.1 ± 0.3 


0.8 ± 0.4 


Ha 




0.05 R 




-1.0 ± 1.5 


4.6 ±1.6 


-1.1 ± 3.3 


0.7 ± 0.2 





9 ± 0.2 


0.2 ± 


.2 


0.4 ± 0.2 


0.3 ± 0.3 





Table 8. Same as Tabled but for correlations with the DIRBE 100 and 240 fim maps. These numbers were obtained from a simultaneous 
fit with the Ha map. Units are fiK, except for the rms of the DIRBE maps {ctdirbe) which is given in MJy sr — . 



DIRBE 


^ DIRBE 


CI 






C2 




C3 




WMAP-K 


WMAP-Ka 


WMAP-Q 


WMAP-V 


WMAP-W 


















\b 


> 30° 










lOO^tm 


0.11 


7.4 ± 


1 


.1 


7.5 ± 1 


.1 


6.5 ± 2.3 




2.9 ± 0.2 


0.5 ± 0.2 


0.0 ±0.1 


-0.4 ± 0.2 


-0.5 ± 0.2 


240^m 


0.27 


6.0 ± 


1 


.1 


3.4 ± 1 


.1 


6.5 ± 2.4 




2.1 ± 0.2 


0.3 ± 0.2 


0.1 ± 0.2 


-0.4 ± 0.2 


-0.4 ± 0.2 


















\b 


> 40° 










100 nm 


0.07 


3.7 ± 


1 


.4 


5.3 ± 1 


.4 


1.8 ± 2.8 




3.7 ± 0.2 


2.4 ± 0.2 


1.4 ± 0.2 


1.4 ± 0.2 


1.3 ± 0.3 


240/J-m 


0.22 


3.8 ± 


1 


.3 


0.7 ± 1 


.4 


-0.8 ± 2.9 




0.9 ± 0.2 


0.2 ± 0.2 


-0.2 ± 0.2 


-0.5 ± 0.2 


-0.3 ± 0.2 


















\b 


> 50° 










100 


0.06 


6.1 ± 


1 


.5 


5.1 ± 1 


.6 


-3.1 ± 3.3 




3.9 ± 0.2 


2.8 ± 0.2 


2.2 ± 0.2 


2.3 ± 0.2 


2.1 ± 0.3 


240/^m 


0.21 


7.1 ± 


1 


.5 


2.0 ± 1 


.6 


-5.3 ± 3.2 




1.3 ± 0.2 


0.6 ± 0.2 


0.5 ± 0.2 


0.1 ± 0.2 


0.2 ± 0.3 



our data and the consistency of our technique to build the 
maps. 

Diffuse galactic synchrotron is not expected to be signif- 
icant over the angular scale to which COSMOSOMAS is sen- 
sitive due to the approximate power law for the galactic sig- 
nal in spherical harmonics of Ci oc l~ 2 . We can use the up- 
per limit on galactic synchrotron at 10 GHz of ATf < 28 [iK 
from the Tenerife results (Gutierrez et al 2000) at high galac- 
tic latitudes to estimate the expected rms in the C2 map of 
< 8 jj,K , assuming a synchrotron spectral index of -3. This is 
in very good agreement with the common signal found with 
the 408 MHz 'cleaned map' (see Tables 7 and 9). As we can 
see in Table 7, the typical amplitudes of the correlated sig- 
nal with the uncleaned 408 MHz map are higher (~ 12 — 17 
/J.K). The excess of correlated signal with respect to the 
desourced map is mostly independent on galactic latitude, 
as it is expected for a random population of extragalactic 
radiosources. 

Comparing the correlated signal of COSMOSOMAS 
with the WMAP.K and WMAP.V channels it is possible to 
estimate the amplitude of the foreground contamination in 
each COSMOSOMAS channel; for instance in the C2 chan- 
nel this foreground signal is 22.5 ± 1.7 /iK at |6| > 30° . 
Galactic synchrotron, free-free and DIRBE correlated emis- 
sion can account for ~ 12/iK. The most likely contributor 



to the remaining foreground signal (~ 19 /iK) is an extra- 
galactic background of radiosources. In fact, although the 
contribution of the strong radio sources has been masked, 
a signal due to extragalactic unresolved radiosources is also 
expected. The amplitude of such contribution can be esti- 
mated from the work by Franceschini et al. (1989), Toffolatti 
et al. (1998) and, more recently, de Zotti et al. (2005). This 
requires the application of the instrument window function 
to the angular power spectrum of the radio sources. Due to 
the Poissonian distribution of sources the power spectrum 
is characterised by a constant Ci. At 15 GHz for sources 
fainter than 0.7-0.9 Jy the value of Ci is expected to be 
in the range 0.75 - 1.25 fiK 2 from the model of de Zotti 
et al. respectively (Toffolatti, private com.). Applying the 
COSMOSOMAS window function calculated in section 2.3, 
we find the range of expected rms in the COSMOSOMAS 
C2 map of 20 - 26 fiK, in agreement with the excess signal 
found. We conclude that apart from CMB fluctuations, the 
signal due to unresolved radiosources is the main contribu- 
tor to the common signal between COSMOSOMAS and the 
lowest frequency WMAP channels. 

We show in Table 7 that the correlation between Ha 
map and the COSMOSOMAS channels is below ~ 5 /xK. 
At 12.7GHz the expected free-free to Ha emission ratio is 
40 /iK/R for an ionized medium with 8000 K (Dickinson et 
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Figure 9. Amplitude of the correlated signal between COSMO- 
SOMAS and WMAP with the DIRBE map at 100 /im for a galac- 
tic cut of |6| > 30°. 



Table 9. Spectral index (in temperature) obtained for three 
galactic cuts from the correlation analysis of the 408 MHz 'cleaned 
map' and the COSMOSOMAS maps. These numbers were derived 
from the coupling constants as a = (i//408 MHz)! 3 . 



Gal. cut 


CI 


C2 


C3 


|6| > 30° 


-3 16+ 003 


-3 n5+°- 04 


, nl +0.08 

J - ul -o.n 


|6| > 40° 


o 9n +0.05 


-3 n3 +0 04 


_o 94+0-07 


\b\ > 50° 


-3 S0+ 07 


-3 16+ 007 
°- lo -0.09 


-3 05+ ' 11 
o.uo_ jo 



al 2004). To estimate the corresponding correlation signal 
for a particular region one multiples this ratio by the rms of 
the template map in that region. For |6| > 30° the rms in 
the Ha map is 0.07 R, which leads to an estimated free-free 
contribution of ~ 3 /iK. This is comparable to the corre- 
lations measured at this galactic cut thus the amplitude of 
the COSMOSOMAS/ WMAP correlated signals with Ha is 
compatible with a free-free spectral index. 

We find clear correlated signals between the COSMO- 
SOMAS channels and the DIRBE maps at 100/im and 
240/im. At \b\ > 30° the weighted- average amplitudes are 
7.3 ± 0.7 fiK and 5.0 ± 0.7 /xK, respectively. The ampli- 
tude of the common signal between COSMOSOMAS and 
DIRBE maps slightly decreases with galactic latitude. The 
cross-correlations of DIRBE with WMAP also shows a pro- 
gressively increasing signal with decreasing frequency with 
some evidence for flattening below ~ 17 GHz (see Figure 9). 
In fact, for the signal correlated with DIRBE 100 fj.ru the 



spectral index between 33 GHz and 22 GHz is -4.3, increases 
to -2.7 between 22 and 16.3 GHz, and becomes nearly flat 
at the COSMOSOMAS frequencies. This spectral behaviour 
is not compatible with the properties expected for thermal 
dust, synchrotron or free- free galactic foregrounds. It is how- 
ever compatible with predictions for anomalous microwave 
emission related to spinning dust (Draine & Lazarian 1998). 
In Watson et al 2005 we find that a signal with similar prop- 
erties is located in a relatively small area (~ 3° x 2°) at 
b — —18.5° located in the Perseus molecular cloud complex. 
To determine whether the DIRBE correlated signal at high 
galactic latitudes found in this paper is associated with simi- 
lar sources or corresponds to a diffuse foreground will require 
better quality data and observations at lower frequencies. 



6 CONCLUSIONS 

The first instrument of the COSMOSOMAS experiment has 
provided temperature maps of the sky emission at frequen- 
cies 12.7, 14.7 and 16.3 GHz covering ~ 9000 square degrees 
with sensitivities from 50 to 120piK per beam. Known point- 
like radio sources with fluxes above ~ 0.8 Jy are detected 
in the three COSMOSOMAS maps. The main results ob- 
tained in the cross correlation analysis of these maps with 
templates of CMB and galactic foregrounds are: 

• The correlated signal between the COSMOSOMAS 
maps and the highest frequency WMAP.V map has an am- 
plitude of 29.7±1.0 fj,K (\b\ > 40°). This is in agreement with 
the expected amplitude of the CMB fluctuations observed 
through the COSMOSOMAS strategy, using the best fit A- 
CDM model to the WMAP power spectrum. 

• The amplitude of the COSMOSOMAS- WMAP corre- 
lated signal increases as one goes to lower WMAP frequen- 
cies, indicating the presence of a foreground signal superim- 
posed over the CMB fluctuations. 

• We find evidence for a contribution of unresolved ra- 
diosources to the COSMOSOMAS maps with an amplitude 
of ~ 19 fiK at 14.7 GHz which we interpret as the result of a 
randomly distributed population of unresolved radiosources. 

• At high galactic latitudes significant correlations for 
COSMOSOMAS and the 408 MHz and 1420 MHz galactic 
templates have been found. The values of the synchrotron 
spectral index are between -3.20 and -2.94 in the 408 MHz- 
16.3 GHz range. 

• The correlation analysis between the COSMOSO- 
MAS/WMAP maps, and DIRBE channels at 100 and 240 
/im indicates also the presence of common signals at high 
galactic latitudes (|6| > 30°). The amplitude of these sig- 
nals rises with a very steep index from ~ 40 GHz up to ~ 20 
GHz, flattening in the range 13-17 GHz with an amplitude 
of 5-7/xK. This behaviour is compatible with the predictions 
of spinning dust models. The small correlation (~ 3 pK) be- 
tween COSMOSOMAS and the Ha map rules out free-free 
emission as the emission mechanism responsible of the dust 
correlated signal detected in these regions. 
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